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Abstract. Numerical simulations are performed of a test scalar field in a spacetime 
undergoing gravitational collapse. The behavior of the scalar field near the singularity 
is examined and implications for generic singularities are discussed. In particular, 
our example is the first confirmation of the BKL conjecture for an asymptotically flat 
spacetime. 



Examining Gravitational Collapse With Test Scalar Fields 



2 



1. Introduction 

Gravitational collapse results in black holes which contain singularities. The 
collapse process is not described by an analytical solution of the field equations, and 
the gravitational field is so strong that perturbation theory is not helpful either. So 
numerical simulation is needed. But even here, the full process of gravitational collapse 
is not usually examined: some simulations examine the formation of black holes but 
not singularities [1] while others examine singularities but not black holes [21 E]. Note, 
however, that there has been some recent progress in numerical simulations capable of 
examining the formation of both black holes and singularities. [5] In the case where black 
holes are examined, the simulation sometimes ends when the black hole forms [1] or the 
simulation goes on, but the interior of the black hole is excised |3]. In the case where 
the singularity is examined, one usually imposes periodic boundary conditions, which 
means that what is studied is singularity formation in a closed cosmology. 

There is a longstanding conjecture due to Belinskii, Khalatnikov, and Lifschitz 
(BKL) as to the general behavior of spacetime singularities. [6] Essentially the BKL 
conjecture is that as a generic singularity is approached in a comoving coordinate system, 
the field equations will be dominated by terms containing derivatives with respect to 
time, and all other terms can be neglected. The BKL conjecture has been verified by 
numerical simulations [21 13] for the case of singularity formation in a closed cosmology. 
Since BKL dynamics are local, it seems likely that the BKL conjecture is also true in 
the case of singularity formation inside a black hole (at least for part of the singularity, 
see [7] for a conjecture about another part of the singularity). However, one would like 
to verify the BKL conjecture by a treatment of the singularities inside black holes. 

One of the earliest studies of gravitational collapse, that of Oppenheimer and Snyder 
[8] was capable of also examining the singularity. In [8] gravitational collapse is modeled 
as a spherically symmetric spacetime containing two parts: the Schwarzschild metric on 
the outside and a Friedmann-Robertson- Walker (FRW) spacetime with non-relativistic 
matter (dust) on the inside. The Schwarzschild and FRW spacetimes have two different 
coordinate systems which must be matched along the boundary. This difficulty can 
be overcome [9] through a coordinate transformation which allows the full spacetime 
to be covered by a single coordinate system. However a more serious difficulty is that 
pressureless dust is not a good model of how matter behaves in the extreme conditions 
of singularity formation. So the properties of dust singularities tell us very little about 
the properties of generic singularities. 

Ideally, one would like to check the BKL conjecture by treating the gravitational 
collapse of a more realistic matter model, like a scalar field, using a single coordinate 
system that follows the formation of the black hole and the singularity. As a preliminary 
step in that direction, we perform numerical simulations of the behavior of a test massive 
scalar field on the spacetime of [8| using the coordinates of [9]. The relevant equations 
and numerical methods are described in section 2, while results are presented and 
discussed in section 3. 
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2. Methods 

A spherically symmetric spacetime can be put in the form [TOl ITT] 

ds"^ = -(1 - ^^)dt'^ + 2ipdtdr + r'^{de'^ + sin^ed^'^) (1) 

Here, r is the usual area coordinate, while t is a time coordinate chosen so that 
V^tVat = —1- We now show how to put the Schwarzschild and FRW metrics in this 
form. The Schwarzschild metric in the usual coordinates is 

ds^ = -(l-'^-M\dt'+(l- —\ ~\r^ + r\de^ + sin^e V) (2) 



We wish to express t as a function of t and r so that the metric takes the form of eqn. 
([1]). Some straightforward but tedious algebra shows that the expression 



t = t-2 



V2Mr + Mln 



yields a metric of the form of eqn. ([1]) with 



(3) 



^ = \/— (4) 
The spatially flat dust FRW spacetime has the metric 

ds^ = -dP + a^{t){dr^ + r^[de^ + sm^Od^^]) (5) 
with the scale factor a given by 

a = KP/^ (6) 

where K is a. constant. We wish to perform a coordinate transformation that gives this 
metric the form of eqn. ([1]). Comparing the angular parts of the metric, it is clear that 
we must have r = ar. It is then straightforward to show that with i = t the metric of 
eqn. ([5]) takes the form of eqn. ([1]) with 

^ - -I m 

The collapse takes place at negative times and is complete at t = 0. The boundary 
between the interior FRW part of the spacetime and the exterior Schwarzschild part is 
given by r = rh{t). for some function rh{t). Since the metric, and therefore ip must be 
continuous across the boundary, we have 



1/3 



n = [IMP) (8) 

Thus the spacetime metric is given by eqn. ([1]) with ijj given by eqn (jT]) for r > r?, and 
ip given by eqn. ([71) for r < r?,. 

We wish to simulate the Klein-Gordon equation on this spacetime assuming 

spherical symmetry for the scalar field. That is, we want a scalar field 0(r, t) that 
satisfies 

VaVy - m^(j) = (9) 
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For a metric of the form of eqn. ([T]) the Klein-Gordon equation for a spherically 
symmetric scalar field becomes 

r 

+ (1 - ilj^)^dr{r^dr4>) - 2iP{drij)dr4> - (f) (10) 

For the numerical method used, we put this equation in first order form: we introduce 
the quantities P and 5* given by 

P = dt<P (11) 
S = 3,4 (12) 

From the definitions of P and S it immediately follows that 

dtci> = P (13) 
dtS = drP (14) 

while eqn. (fTOl) becomes 

2th 

dtP = 2i)drP + {dti))S + {dri))P + —P 

r 

+ (1 - iP^)^dr{r^S) - 2i){dri))S - m^(P (15) 

The equations to be evolved are eqns. f lT3lfT5|) . We evolve using finite differences, 
where each scalar /(r) is represented by the quantities fi = f{iAr). Spatial derivatives 
are evaluated using standard centered differences, that is 

drf ^ ~ (16) 
2Ar ^ ' 

However, to maintain accuracy near the origin, the quantity r"^9r(r^S') requires special 

treatment. We define 5"+ and by 5+ = (S'j+i + Si)/2 and = {Si + S'j_i)/2. 

Correspondingly we define r_|_ = rj + (Ar/2) and r_ = — (Ar/2). Then we use the 

finite difference 

l9.(r^^)^-^(r^5^-r^5_) (17) 

The quantity ijj and its derivatives are evaluated analytically using the formulas in eqns. 
01]) and ([7]). Eqns. ( ITBlfT^ are evolved using a three step iterated Crank- Nicholson 
method [12] with Kreiss-Oliger dissipation [T3] . 

At the origin, (gridpoint i = 1) smoothness requires that 5* vanish and that and 
P have vanishing derivative with respect to r. We implement this condition as 

01 = (402 - 03)/3 (18) 

Si = (19) 
Pi = (4P2 - P3)/3 (20) 
At the outer boundary, we impose the condition that 0, S and P all vanish. 
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3. Results 

In the simulations, let be the number of gridpoints and Tmax be the maximum value 
of r. Then we have Ar = r^g,^/{N — 1). = 800 was used for the runs resulting in 
figures [lUHl = 16000 was used for the runs resulting in figures \MTli We choose initial 
data for the scalar field at a time to to take the form P = and 

with S given by the derivative of this expression. Here, tq and a are constants. For all 
the results displayed, the parameters are to = —20, rmax = 80, ro = 35, a = 1.5, and 
M = 5 where M is the mass of the Schwarzschild black hole. This initial configuration 
of the scalar field is shown in figure [H We run the simulations to a time near the 
singularity. Figures [2] and [3] show the configuration of the scalar field near this final 
time. Both figures are for t = —0.23 with m = for figure |2]and m = 1 for figure |3l 
The initial scalar field can be thought of as a linear combination of a left mover and a 
right mover. The right mover escapes from the black hole, while the left mover plunges 
towards it. Some part of the left mover scatters off of the black hole, while another part 
is captured and approaches the singularity. 

We are particularly interested in finding the behavior of the scalar field as the 
singularity is approached. In our coordinate system, the singularity is approached as 
t — )■ but only in the FRW portion of the spacetime. To examine the scalar field 
behavior near the singularity, we plot the scalar field at three different times near t = 
and deduce a trend from those plots. The scalar field (j) is plotted as a function of r 
in a small region near r = for the times —0.02, —0.01 and —0.005 for the m = 
case in figure H] and for the m = 1 case in figure [5l It is clear that the trend is that 
as the singularity is approached, the scalar field gets a larger amplitude and a steeper 
profile. To see whether this behavior can be understood quantitatively, we turn to the 
BKL conjecture: the claim that as a generic singularity is approached in a comoving 
coordinate system, the field equations will be dominated by terms containing derivatives 
with respect to time, and all other terms can be neglected. To apply this conjecture, 
we consider the Klein-Gordon equation (Q in the untransformed FRW metric and 
coordinates given by eqns. ([5]) and ([6]). Keeping only terms containing time derivatives, 
we find 

i-%{Pdi^) = 0. (22) 
The general solution of this equation is 

= r^A{r) + B{r) (23) 

for some functions A{r) and B{r). In particular, the BKL conjecture leads us to 
expect that as the singularity is approached, t(j) should approach a function of f. Note, 
that expressed in terms of the t, r coordinates of our simulation, we have t = t and 
f = K^^t^'^^^r. Thus, we are led to expect that as the singularity is approached, if 
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Figure 4. ^ r for m = at i = -0.02 (solid), -0.01 (dashed) and -0.005 (dotted) 
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Figure 5. vs r iovm = 1 at t = -0.02 (solid), -0.01 (dashed) and -0.005 (dotted) 




Figure 6. t(j) vs t ^/^r for m = at < = 
(dotted) 



-0.02 (solid), -0.01 (dashed) and -0.005 
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Figure 7. tcj) vs f'^/^r for m = 1 at i = -0.02 (solid), -0.01 (dashed) and -0.005 
(dotted) 
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Figure 8. (f> vs r ior m = at t = -0.0055 (solid), -0.0030 (dashed) and -0.0017 
(dotted) 

we plot t(f) vs f'^^^r the plots should approach some limiting profile as t — )■ 0. To see 
whether that expectation is realized, such plots are given in figure |6] for m = and [7] 
for m = 1. These figures are consistent with t0 approaching a limit as the singularity is 
approached. However, in order to be more definite, we need to approach the singularity 
more closely. The results of simulations which show the scalar field at times closer to 
the singularity are shown in figures [HlfTTl Here figure [8] shows cf) vs r for m = and the 
times —0.0055, —0.0030, and —0.0017. Figure [9] shows ws r at the same times for the 
m = 1 case. The corresponding rescaled quantities, that is tcp vs t~'^^^r are shown in 
figure [To] for m = and in figure [TT] for m = 1. 

It is clear from figures [TO] and [TT] that tip approaches a limit as the singularity is 
approached. Since this is just what one would expect from the BKL conjecture, our 
simulations lend support to the notion that the BKL conjecture holds, not only in the 
extensively studied case of closed cosmologies, but also in the more physically interesting 
case of asymptotically fiat spacetimes. To explore this issue further, it would be helpful 
to go beyond the case of a test scalar field on a background spacetime. In particular. 
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Figure 9. (p vs r iov m = 1 at t = -0.0055 (solid), -0.0030 (dashed) and -0.0017 
(dotted) 
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Figure 10. t(j) vs t-'^'^r for to = at t = -0.0055 (solid), -0.0030 (dashed) and 
-0.0017 (dotted) 
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Figure 11. vs t'^/^r for to = 1 at t = -0.0055 (solid), -0.0030 (dashed) and 
-0.0017 (dotted) 
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one could treat the collapse of a self-gravitating scalar field using the method of[S] and 
see whether the singularity produced in the collapse satisfies the BKL conjecture. 
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